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Abstract: The Passive Advanced Unit (PAU) for ocean monitoring is a new type of 
instrument that combines in a single receiver and without time multiplexing, a polarimetric 
pseudo-correlation microwave radiometer at L-band (PAU-RAD) and a GPS reflectometer 
(PAU-GNSS/R). These instruments in conjunction with an infra-red radiometer (PAU-IR) 
will respectively provide the sea surface temperature and the sea state information needed 
to accurately retrieve the sea surface salinity from the radiometric measurements. PAU will 
consist of an array of 4x4 receivers performing digital beamforming and polarization 
synthesis both for PAU-RAD and PAU-GNSS/R. A concept demonstrator of the PAU 
instrument with only one receiver has been implemented (PAU-One Receiver or PAU-OR). 
PAU-OR has been used to test and tune the calibration algorithms that will be applied to 
PAU. This work describes in detail PAU-OR' s radiometer calibration algorithms and their 
performance. 

Keywords: radiometer, reflectometer, pseudo-correlation, calibration, field-programmable 
gate array (FPGA). 



1. Introduction 

Sea surface salinity can be indirectly measured through the variations of the brightness temperature 
due to the change of the sea water dielectric constant with respect to temperature and salinity [1]. 
However, the brightness temperature also depends on the sea surface roughness. Despite the field 
experiments performed in the past years to improve our understanding of this effect [2,3,4], the sea 
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surface roughness correction still remains one of the most critical corrections needed to retrieve the 
salinity with the required accuracy since it cannot be simply parameterized in terms of the wind speed 
and/or the significant wave height. Inaccuracies in the brightness temperature direct model may induce 
significant errors (biases) in the retrieved salinity [5]. 

The Passive Advanced Unit (PAU) (Fig. 1) is a new concept instrument originally proposed in 2003 
[6]. PAU consists of a suite of three instruments operating synergistically: 

• PAU-RAD: a new type of pseudo-correlation L-band radiometer to measure the brightness 
temperature of the sea surface, 

• PAU-GNSS-R: a GPS-reflectometer to measure the sea state, and 

• PAU-IR: an IR radiometer to measure the sea surface temperature. 

Figure 1. PAU overview. 




PAU-RAD is a full polarimetric radiometer, so that the delivered results are the four Stokes 
parameters (eqn. 1). The Stokes parameters fully characterize the polarization state of partially 
polarized thermal radiation. This polarized thermal radiation is related to different phenomena such as 
the angle incidence, the sea surface salinity (SSS), the sea surface temperature (SST), the sea surface 
roughness (SSR), wind direction, foam... 
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where: 



• I,Q,U and V are the four Stokes parameters, and 

• E v and E h are the electric fields at vertical and horizontal polarizations. 
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The first and the second parameters (7, Q) are the sum and the difference of the brightness 
temperatures at vertical and horizontal polarizations. It is known that the sea surface salinity can be 
retrieved from I and Q [1]. The third and the fourth parameters (U, V) characterize the correlation 
between these two orthogonal polarizations. Wind direction information can be retrieved from U and V 
[7]. 

The GPS-reflectometer uses the same RF-IF front-end as the L-band radiometer, although the digital 

beamformer synthesizes beams pointing towards the specular reflection points of the GPS signals [8,9]. 

Ar 

The relationship of the standard deviation of the SSR and the brightness temperature is obtained ( — - ) 

&SSR 

[10] from this instrument. Finally, the third element is a commercial 8-14 um thermal IR radiometer 

[11], from which the SST is obtained. PAU [12] will retrieve the SSS from radiometer outputs (Stokes 

Ar 

parameters) and correcting the sea state effect using information from the reflectometer ( — -) and 

&SSR 

from the IR-radiometer (SST). 

PAU-RAD consists of an array of 4x4 dual polarization receivers integrated behind a patch antenna 
[13], whose outputs are digitized, and then properly calibrated and combined to produce several beams 
using a digital beamformer [14]. In this implementation the PAU-RAD operates at the GPS LI 
frequency, sharing the RF-IF front-end with the PAU-GNSS/R. PAU-RAD' s digital design has been 
split in two Field Programmable Gate Array (FPGA) [14]: one is called Control Unit (CU) and 
manages the system and the another one is called Arithmetic-Logic Unit (ALU) and implements the 
radiometer algorithms. 

Section 2 presents the receiver design and operation. Section 3 describes the calibration algorithms 
implemented, which includes the injection of correlated noise at two different levels, and a Kalman 
filter [15] to better track the phase drifts. Calibration results are presented at the end of this section. 
Finally, section 4 presents the main conclusions and future research lines. 

2. PAU-RAD design overview 

Each element in the array has a dual-polarization antenna (horizontal and vertical) and each 
polarization is thereafter divided in two using a Wilkinson power splitter. The PAU radiometer idea is 
that the output voltage signals are proportional to the antenna signal plus or minus the noise added by 
the resistor of the Wilkinson power splitter. A detailed noise wave analysis of this novel topology is 
described in [12]. The receiver topology is described in detail in [13]: The input signals are 
demodulated at an intermediate frequency (IF) of 4.309 MHz and have a 2.2 MHz bandwidth. The 
analog signals are then digitalized at 8 bits at a sampling frequency of 5.745 MHz. 

To avoid the mechanical scan of the antenna, a digital beamformer with the 4x4 element array has 
been designed so that the beam can be steered up to +20° from the array boresight (45° incidence angle) 
in 5° steps (incidence angle from 25° to 65°). The PAU-RAD basic requirements are summarizes 
below: 
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• I/Q demodulation of 64 input signals (16 elements x 2 polarization/element x 2 receiving 
chains/polarization), 

• maximum amplitude unbalance within +1 dB. Residual amplitude unbalance better than 
+0.1 dB. Limited by the number of bits, 

• phase calibration between receiving chains better than 2° rms (random initial error within 0° 
and 360°), 

• computation the four Stokes parameters using correlation techniques, 

• system re-configurability, so that the integration and calibration times, and the inter- 
calibration period may be optimized as a function of the scenario and application, and 

• simultaneous generation of two simultaneous beams, with the computing resources 
available. 

As it has been said in the introduction, the reflectometer and the radiometer share the same RF-IF 
front-end, a commercial GPS down-converter. As a result both, radiometer and reflectometer, operate 
at 1575.42 MHz. However, it is known that for the best performance, radiometer should operate at the 
1400-1427 MHz band and reflectometer at 1575.42 MHz. This non-optimal scenario is assumed 
because the aim of this prototype is to demonstrate the PAU concept and the scenario allows reducing 
the hardware resources required. On the other hand, we have to deal with some handicaps. The fist is 
the radiometer uses the 2.2 MHz GPS bandwidth, instead of the about 20 MHz available in L-band. 
Obviously, it has an impact on the radiometric sensitivity that can be compensated by increasing the 
integration time. The second handicap is the contamination that GPS signal can introduce in the 
radiometric measurements, which has to be quantified despite of the spread spectrum modulation GPS 
signal is lower than the noise level. Total GPS power density received over the Earth's surface is about 
-153 dBW/m z . The power received after a reflection over the Earth's surface depends on the receiver's 
height and on the antenna' s effective area. In the worst case, where the specular reflection point is at 
the antenna's boresight and considering a surface reflection coefficient of T=0.7, a height of h=100 m 
and an antenna's effective area of A ef ^0.16 m the radiometric measurement contribution is AT= 
P r //c B B=2.3 K (~ 4.6 psu salinity error in warm water). However, if the reflection comes through the 
array pattern with an attenuation at least of 15 dB, the radiometric measurement contribution is just 
AT= 0.03 K (~ 0.06 psu salinity error in warm water), below the 0.1 psu requirement of GODAE for 
open ocean [16]. Taking into account that the radiometer beam can be electronically steered, the 
chances of having a significant interference level are negligible. 

PAU-RAD is a proof of concept of the PAU-SEOSAT/INGENIO [17] instrument that has been 
selected for the B phase of SEOSAT/INGENIO (stands for Spanish Earth Observation Satellite, 
Spanish contribution to GMES). satellite. In this PAU version the radiometer will operate in the 1404- 
1423 MHz band (bandwidth of 19 MHz) and the reflectometer at the LI GPS frequency, using 
common RF chains and two separated IF chains. 
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2. 1. PAU-RAD analog receiver overview 

Each receiver (Fig. 2) has two linear polarization antennas and two chains per polarization (four chains per 
receiver). 

Figure 2a. PAU's RF part of the Figure 2b. PAU's IF part of the 
elemental receiver with box. elemental receiver with box. 





A switch follows the antenna and commutes between the antenna signal (SA) and the correlated noise 
(SCN) (Fig. 3) from a controlled noise source (CNS), used for calibration purposes. A CNS has been 
designed to provide three temperature points (300, 400 and 500 K), so as to perform two types of 
calibration: internal hardware calibration and absolute radiometric calibration. 

Figure 3. Analog receiver configuration for hardware calibration by noise injection. 
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After the switch, the Wilkinson power splitter divides the signal in two, one for each receiver chain. The RF 
stage amplifies the signal by 33 dB. At the downconverter output the signals are centered at an IF of 4.309 
MHz with a bandwidth of 2.2 MHz. The last stage is formed by the adjustable gain video amplifiers and an 
impedance matching network. From here the signals are sampled using 8 bits analog-digital converters 
(ADC). The fine amplitude adjustment is then performed in the digital part, using advanced processing 
algorithms as described in section 3. 
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2.2. Analog Acquisition 

The number of bits used is a trade-off between the radiometric sensitivity of the full array, which 
depends on the clipping effect [18] and on the effective number of bits, and FPGA resources. The 
analysis developed on [19] shows that the best performance is given when the ADC's conversion 
window and the standard deviation of the signal is Vpp/csiGNAL = 9.09, where V pp means the window's 
voltage peak to peak. The same work analyzes the performance as a function of the number of bits, the 
result is that 8 bits is enough to neglect the digitization error and to ensure a main beam efficiency 
larger than 96 % of the digital beamformer. 

Considering these results, the analog signals are sampled at 8 bits and 5.745 MHz. The sampling 
frequency is higher than the IF which allows a kind of "band-pass sampling" [20] downconverting the 
input signals at 0.25 times the digital frequency (1.436 MHz). Digitization of the 4 receiver chains (4 
per element) is performed with 2 dual-channel 8-bit ADCs, which share the output bus. An internal 
multiplexer determines which digitalized channel is present at the ADC outputs. It is controlled by an 
1 1.49 MHz signal, so in half -period of CLK each ADC is able to put its channels in the output bus. The 
control signals are all supplied by the FPGA to keep the synchronism with the internal digital 
processes. 

Moreover, according to the relationship V pp = 9.09- csiGNALthe signals amplitude have been adjusted 
to 110 mVrmsto clip less than 0.1% samples when the ADC dynamic range is set to 1 V. The digital 
output values go from 0 (-0.5V) to 255 (+ 0.5V), so a 128 offset value is present at the output data, that 
will have to be subtracted afterwards. 

The signals from each element enter into the FPGA multiplexed in time by a factor of 4 (4 times 
faster). Since there are 16 elements, the processing unit must process 367.68 Mbytes of data per 
second. To process this data rate properly the hardware must operate, at least, at 4 times the sampling 
frequency (22.98 MHz). 

2.3. PAU-RAD digital overview 

The digital part of PAU-RAD has been split in a Control Unit (CU) that manages the system and 
another block called Arithmetic-Logic Unit (ALU), where the radiometer algorithms are implemented. 

2.3.1. PAU-RAD Arithmetic Logic Unit (ALU) 

The PAU-RAD ALU processes the acquired data and finally computes the four Stokes parameters 
[21]. The digital implementation to compute the first two elements has been performed as follows: 




(2a) 




(2b) 
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where: 



• S\ and S2 are the signals from the vertical antenna with the Wilkinson noise (Wilkinson noise 
is defined as the noise generated in the Wilkinson load) and the receiver noise, (Si ,2 signals 
are deeply analyzed on 3.4), 

• S3 and S4 are the signals from the horizontal antenna with the Wilkinson noise and the 
receiver noise, (Si^ signals are deeply analyzed on 3.4), 

• A,2,3,4 and <2i,2,3,4 are the in-phase and quadrature signals for Si ; 2,3,4 signals, respectively, and 

• Tj is the integration time. 

And the digital implementation to compute the second two elements has been performed as 
follows: 



• Si and S2 are the signals from the vertical antenna with the Wilkinson noise and the receiver 
noise, 

• S3 and S4 are the signals from the horizontal antenna with the Wilkinson noise and the 
receiver noise, 

• 7i ;3 and Q l 3 are the in-phase and quadrature signals for Sjj signals, respectively, and 

• T i is the integration time. 

In order to save FPGA resources and maximize their use, signal time-multiplexing and hardware 
reuse techniques have been intensively used. This increases the synchronism complexity because it is 
necessary to relate the ADC output to a synchronous signal to know which signals are being processed 
at each clock cycle. Furthermore, to improve the maximum operating frequency two techniques are 
used: 

• parallelization, that involves the use of sequential circuits (changes a multiple cycle serial 
operation to a multiple stage parallel operation), and 

• pipelining, that involves the use of combinational circuits (changes a high delay 
combinational circuit to multiple low delay circuits with registers). 

In these conditions the digital design operating frequency must be the same as the sampling 
frequency. Actually, data is processed at twice the sampling frequency in order to obtain a second 
steered beam. 




(2c) 




(2d) 



where: 
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The data path is schematically shown in Fig. 4. A correction must be first applied to the input 
signals to adjust their level before being processed. This is related to the output of the ADC that gives a 
binary quantification of the analog signal between 0 and 255. It is therefore necessary to subtract 128 to 
all the signals while, at the same time, converts the signals into 2's complement. 

Figure 4. PAU-ALU-RAD diagram. 
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In order to perform the beamforming and the pseudo-correlation, all channels must have a zero 
relative phase and the same amplitude. Phase errors are largely related to the fact that the 64 analog 
receivers do not share the same LO. A 10 MHz clock reference is distributed to all receivers' PLLs, 
which lock with a random phase what is uniformly distributed within 0 to 2n radians. Amplitude 
errors are introduced by gain variations (aging and drifts) from each receiver's chain, and are adjusted 
within +1 dB. 

As discussed in section 2.2, the incoming signals are centered at 0.25 times the digital frequency. 
Then, it is possible to simultaneously down-convert, apply the beam steering phase, and perform the 
phase and amplitude calibration of the signals with a single multiplication between the signal and the 
digital local oscillator signal: 

S L0 (n) = aexp(-j(w n n-0 beam ))with n=l,2,3,... ( 3 ) 



where: 



w n = 0.25 times the digital frequency, 

a : complex calibration coefficient (amplitude and phase) and 
e beam : beam steering phase. 
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To synthesize this in the digital domain, a Multiplexed Digital Down-Converter (MDDC) block has 
been developed (Fig. 5) which is also fully adaptive to create new beams or dynamically compensate 
the phase/amplitude error correction. The MDDC consists of a phase and amplitude variable Local 
Oscillator (LO), whose values are determined by the PAU-RAD control unit (CU) as a function of the 
beam steering, and the estimated phase and amplitude errors. The LO is updated when a new beam is 
synthesized or the instrument recalibrated. Note that there are 16 MDDCs and 64 LOs, one for each 
receiver. 



Figure 5. MDDC sketch. 
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After the MDDC the digital beamformer signal is formed by adding the 16 signals that came from 
the receivers using a weighted adder. The next step is to low-pass filter the I/Q composed beams 
signals, to do this an ITR filter has been chosen instead of FIR one, due to the smaller order required to 
obtain the same attenuation. Once the signals are filtered, they are input to the pseudo-correlator block 
that multiplies them by pairs (Eqn. 2), the results are accumulated in the integrator block. 

PAU-RAD CU switches the input analog receiver signal from antenna signal to a correlated noise 
signal to estimate the phase and amplitude errors. The correlated noise is needed to determine the 
relative amplitude and phase errors between channels and calibrate them. 

2.3.2. PAU-RAD Control Unit (CU) 

PAU-RAD CU is built on an embedded processor which is implemented in the FPGA. It is a 32 bit 
processor with 16 bit instructions specially developed to be programmed in C language. PAU-RAD CU 
has three main functions: 

• it is both, the interface with the PAU-RAD ALU and with the external interface of PAU- 
RAD; controlling the system operation mode and transmitting to PAU-RAD ALU the state, 

• computes the corrections to be applied at the digital LO to calibrate the signals, and 

• attends the remote control based on a TCP/IP connection. 
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PAU-RAD CU has been designed as a finite state machine, which is always working in one of the 
different states according to the radiometer' s operation mode at each time instant. The block diagram in 
Fig. 6 shows the different states of PAU-RAD CU and the transitions between them. There is a C main 
function to implement this State Diagram, and a subroutine to each state, therefore the descriptions of the 
different states are an explanation of the corresponding subroutine features. There are two types of states: the 
extended ones that perform test functionalities (PHASE SWEEP (PS), CLOSE LOOP (CL) and OPEN 
LOOP (OL)), and the fundamental ones that perform the core of the radiometer functions 
(INITIALIZATION (INIT), CALIBRATION (CAL), CORRELATED INTERNAL REFERENCE 
(CIR), UNCORRELATED INTERNAL REFERENCE (UIR) and ANTENNA MEASUREMENT 
(AM)). 

Figure 6. PAU-RAD CU extended datagram. 
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The different states are described below: 

• INIT: initialization state. It has two main functions: 

a. sets the default parameters of the system when resetting and 

b. it listens the remote commands coming from the external interface. In a reset process the 
PAU-RAD CU transmits to the PAU-RAD ALU the system timer. 

• CAL: hardware calibration state. It switches the input signal to the correlated noise source 
(CNS) and starts the phase and amplitude error estimation. The LO is modified using the 
correction values once the error is achieved. 

• CIR: internal radiometric calibration state. It switches the input signal to two well known 
temperature points of the CNS. The signals are correlated and sent to the external interface 
through the TCP/IP connection. These results are used to set a unique relationship between 
the detected power and the brightness temperature. The aim of this state is to reduce the time 
between external hot-cold calibrations. 
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• UIR: offset estimation state. It switches the input signal of each chain to a matched load and 
correlates the signals. The result is the offset correction to be applied at the radiometric 
measures. 

• AM: Computes the Stokes parameters. It switches the input signal to the antennas. It is the 
state in which the system spends most of the time. 

• OL: temporal phase and amplitude drift analysis state, only used for test purposes. It 
switches the input signal to a CNS' predefined point. LOs are set to the same amplitude and 
phase, them the power and the phase of each chain is computed. 

• CL: phase and amplitude calibration performance state, used only for test purposes. It 
switches the input signal to a CNS' predefined point. The system is always estimating and 
correcting the phase and amplitude errors. 

• PS: LO phase sweep analysis state, only used for test purposes. It switches the input signal to 
a CNS' predefined temperature point. The system changes constantly the LO's phase and 
computes the correlation between channels. 

The standard sequence (CAL-CIR-UIR-AM) is only modified if a remote command has been 
received. In this case, the system moves to INIT state. The time spent in each state depends on the 
hardware stability. The variable v_timecal is a timer that can be configured depending on the scenario 
and determines the inter-calibration period. It is the time that the system spends in the AM. Each 
fundamental cycles, except AM, are executed some number of times before the next state. 

Finally, the CU has a remote reset watchdog (the reset signal comes from TCP/IP connection) to 
periodically check the status of the connection, and reset the CU if necessary. 

3. Calibration algorithms implementation and tests 

In order to test and improve the calibration algorithms for PAU a one-receiver version has been 
developed (PAU-OR) [22]. It includes PAU-RAD and PAU-GNSS/R, but not PAU-IR. 

Figure 7. PAU-OR general view. 
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PAU-OR (Fig. 7) has a full polarization seven-patch array antenna with 25° beamwidth and 17.5 dBi 
directivity. Analog circuitry is located inside a thermal isolated and regulated box using Peltier cells 
(achieving a thermal stability, standard deviation of the internal temperature, lower than 0.3°C). A left 
hand circularly polarized signal is composed and it is sent to PAU-GNSS/R. A 15 dB ENR noise 
source is used for the CNS to perform the hardware calibration and for internal reference radiometric 
calibration. 

In PAU-OR, and by extension in PAU-RAD, there are two types of calibration. The first one is 
related to the internal relative calibration between the analog receivers that have different gains and 
phases. Needed to perform the beamforming and operates as a pseudo-correlation radiometer. The 
second one is the radiometric calibration to set a unique relationship between the cross-correlation 
outputs and the Stokes elements. 

The CNS has been designed to achieve these two types of calibration, allowing the generation of 
three different temperature points. 



3.1. Internal relative phase and amplitude calibration 

The analog receiver (Fig. 3) has two possible inputs: one from the antenna (MA mode) and the other 
from the CNS, for calibration purposes (CIR mode). By commuting the input switch to the correlated 
noise input the same signal is applied at the input of each chain, which allows estimating the amplitude 
and phase differences among them. The equations obtained from PAU are extensively developed in 
[12] where all non-ideal terms such as the Wilkinson power splitter unbalance or the imperfect 
matching and isolation are taken into account. Note that in this paper the system is assumed to be ideal 
to simplify the mathematics and to clarify the calibration algorithm. If the system imperfections are 
accounted for, additional terms appear, but cancel out in the differential algorithm described below. 

When the input of the analog receivers is commuted to the CNS path, the four ideal output signals 
are given by: 
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where signals are defined as effective voltage waves: 



• Scn : is the noise signal from the CNS (same for the four channels), 
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• S w i: is the thermal noise signal generated at the first Wilkinson power splitter that divides 
the noise from the CNS in two, one for the upper branch (chains 1 and 2, Fig. 3) and the 
other for the lower branch (chains 3 and 4, Fig. 3), 

• S W 2 and S W 3: are the thermal noise signals from the second and third Wilkinson power 
splitter, at the upper branch and at the lower branch respectively. Note that these signals are 
uncorrelated among them, 

• Sreci, Srec2, Srec3 and S RE c4: are the thermal noise signals introduced by from each 
receiver. Note that these signals are also uncorrelated among them, and 

• a, f3, y: are the complex coefficients that contain the relative phase and amplitude error 

with respect to the first channel. 



To estimate the complex coefficients a, /?, y, the following correlations are obtained using the 
signals defined in Eqn. 4: 



T +T T 
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(5a) 
(5b) 
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(5c) 
(5d) 



where: 



• T C n: noise temperature from the CNS, since T CN □ l^S CN \ j (Note that the equivalence to 

brightness temperature has been normalized by ke). 

• T w i: noise temperature introduced by first Wilkinson power splitter, note that T wl □ T phl □ T ph2 , 

that are the receivers physical temperature and 

• T W 2: noise temperature introduced by the upper branch Wilkinson power splitter, note 
thatr W2 □ T ph2 □ T phl , that are the receivers physical temperature. 



In order to compare Eqns. 5a and 5b with 5c and 5d the ± and 7^, terms that must be 

cancelled. To eliminate these terms a differential calibration approach is used, by injecting two 
different noise temperatures (Tcni and Tcn2), and subtracting the results for each channel: 
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The complex coefficients a, P and y are them obtained without the "a priori" knowledge of the 
receiver's noise temperature. The digital LO is then modified and, as a result, the channels are 
calibrated. 

Note that: 

• Equation 5 and Eqn. 6 the k B terms have been deliberately omitted so that "powers" are 
directly expressed as temperatures (Kelvin). This analysis can be directly extended to the 
case of non-ideal power splitters and receivers as detailed in [12], and 

• As is said, the previous analysis has assumed that the Wilkinson power splitters were ideal, 
and that all inputs were perfectly matched. A detailed analysis was performed in [12] which 
shows that when the differential technique is used, the non-ideal terms are cancelled out, 
thus reducing to Eqn. 6. 

3.2. Test of the internal relative phase and amplitude calibration algorithms 

A set of test has been performed in order to evaluate the algorithms performance. All these tests 
have been done in CL mode, the system is always injecting two noise levels to estimate and correct the 
phase and amplitude errors on a snap-shot basis every 0.53 s. (integration time). Figure 8 shows the test 
set-up. Note that the receiver is outside the thermal isolated box to show the whole test set-up. 

Figure 8. PAU-OR general view, set up for test purposes. 



FPGA PAU-GNSS/R 



Noise source 




Sensors 2008, 8 



4406 



Figure 9 shows the different output power ( in counts: arbitrary measurement units, AMU), during 
the warm up transient, for the four receiving channels. After a warm-up transient of 50 min., the 
amplitude values stabilize, and the calibration algorithm can then be applied. 

Figure 9. One receiver amplitude temporal drifts, without any correction. 
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After applying the calibration algorithm at each measurement. The means and the standard 
deviations for each channel are 011,2,3,4= 58.63 dBAMU, within the goal of +0.1 dB residual amplitude 
calibration error, and ai= 0.017, a 2 = 0.032 a 3 = 0.016 and a 4 = 0.015 dBAMU, respectively. (Fig. 10b) 



Figure 10a. One receiver amplitude 
temporal drifts, without any correction. 



Figure 10b. One receiver amplitude after the 
differential amplitude calibration algorithm. 
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Note that in Fig. 10a the powers are expressed in arbitrary measurement units (AMU, 1 AMU = 1 
count). Before calibration the mean and the standard deviation for each channel are mi= 58.63, m 2 = 
58.53, m 3 = 58.64 and m 4 = 58.15 dBAMU and ai= 0.015, c 2 = 0.016 a 3 = 0.003 and c 4 = 0.003 dBAMU, 
respectively. 

Figure 1 1 shows the phase drift of the four channels just after turning on the receiver. Again, after a 
warm-up transient of 50 min., the phase values stabilize, and the calibration algorithm can be applied. 
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Figure 11. Phase temporal drifts of the four channels, without any correction. 
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Figure 12 shows the residual phase error as the result of the application of the described algorithm 
applied. After the transient period, the residual phase error standard deviations are: oi= 0.069°, 02= 
1.343°, 0-3= 1.116° and c 4 = 1.302°. 



Figure 12a. Residual phase error, after 
the phase calibration. 




Figure 12b. Residual phase error histogram 
after the phase calibration. 
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In order to reduce even further the residual phase error a closed-loop Kalman filter [15] has been 
implemented (Fig. 13). A Kalman filter is an efficient recursive filter that estimates the state of a 
dynamic system from a series of incomplete and noisy measurements, in this case the filter's input is 
the estimated phase error that contains measurement noise. The filter's output is the estimated phase 
error with smaller noise level. 
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Figure 13. Digital hardware calibration for phase error estimation using a Kalman filter. 
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There are two main parameters in a Kalman filter to be tuned: the standard deviation of the noisy 
signal and the gain coefficient. The first one, the measurement of the standard deviation of the noisy 
signal, was experimentally determined through several iterations. The second one, the filter gain, is a 
trade-off between the convergence speed and the residual error and the best results have been obtained 
with a 0.95 gain. 

Figure 14a shows the temporal evolution of the residual phase error using a closed-loop Kalman 
filter. After a transient of about 100 s., when there is some ringing, the residual phase error stabilizes to 
smaller values. 

Figure 14b shows the phase error histograms of the residual phase error using a closed-loop Kalman 
filter. The standard deviations of the phase of each channel are: o"i= 0.012°, 02= 0.227°, 03= 0.278° and 
04= 0.189°, while the means of each of the four channels is zero. These residual standard deviation 
values are very low and allow performing the pseudo-correlation radiometer and the digital beam- 
forming accurately. 



Figure 14a. Residual phase error, using a 
closed-loop Kalman filter. 
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Figure 14b. Phase error residual histogram, 
using a closed-loop Kalman filter. 

Pr:b;D I Ty J sii:ini:ri F.rcion Ciiarnei ■ Probsbihtj j siri DLition Fundcn Charnel 2 




-15 -10 -5 0 5 10 -3 -6 -4 -2 0 2 4 



Sensors 2008, 8 



4409 



Note that the histograms (Fig. 14b) look like those of a sinusoidal histogram due to the sinusoidal 
oscillations of the Kalman transient ringing. 

3.3. PAU hardware calibration global algorithm 

This algorithm is implemented in the PAU-RAD-UC, the algorithm has six steps that can be 
summarized as follows: 

• Switching the input to the CNS input, using the highest temperature point available, and 
acquisition of one set of data from the PAU-RAD-ALU. 

• Data acquisition from the PAU-RAD-ALU in order to estimate the relative phase error. For 
each channel compute and set the new error phase values for the LO. 

• Data acquisition from the PAU-RAD-ALU in order to estimate the amplitude error, first 
acquisition. 

• Modification of the CNS to obtain the medium noise temperature point and acquire data 
from the PAU-RAD-ALU in order to estimate the amplitude error, second acquisition. 

• Compute and set the final values for the LO, taking into account phase and amplitude errors, 
and steering phase. 

3.4. Radiometric absolute calibration using internal references 

The classical method to calibrate a radiometer consists of using (at least) an external T H ot, (eg. 
microwave absorber) and an external T C old, (eg. the clear sky). Using these well-known reference 
points, an unique relationship between the antenna temperature and the measured power is obtained. 
This method requires a mechanical movement of the antenna. To avoid this inconvenience, a 
radiometric calibration based on internal temperature references has been implemented. The CNS 
system generates a third well-known temperature point at the input of the receiver chain. The input 
noise is the same as that of a microwave absorber placed in front of the antenna, both at the same 
physical temperature than the matched load. With these results and using a previous measurement of 
Tsky (which obviously includes the antenna Ohmic losses, r\a) the correction coefficients can be 
obtained. 

4. Conclusions and future work 

The calibration algorithm of PAU-RAD has been presented on a one-receiver demonstrator: PAU- 
OR. The instrument has been implemented, the calibration algorithms tested and the result presented 
and discussed. On one hand, to take advantage of using the same band for both reflectometer and the 
radiometer, its beam must be steered to a direction where the GNSS -reflected signal is attenuated by at 
least 15 dB (AT= 0.03 K). On the other hand, the reflectometer beam must be pointed to the specular 
reflection. 

The digital design can be easily modified to generate two simultaneous beams in real-time. The 
whole design, analog and digital parts, has been satisfactorily tested. Moreover, the PAU-RAD CU has 
been designed in C language over a microprocessor giving it versatility. The state diagram has been 
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developed to obtain the best performance in any kind of scenario and allows changing the main 
radiometer' s parameters. 

The main results are summarized below: 

• Differential amplitude and phase calibration algorithms have been developed in order to 
mitigate the offset introduced by the subsystems non-idealities and to give robustness to the 
error estimation. 

• The results obtained with de differential amplitude calibration algorithm are better than the 
required O.ldB: d= 0.017, c 2 = 0.032 c 3 = 0.016 and c 4 = 0.015 dBAMU. 

• The results obtained with the differential phase calibration algorithm are smaller than 1.5°: 
d= 0.069°, c 2 = 1.343°, c 3 = 1.116° and c 4 = 1.302°, and 

• the improvement over the standard phase calibration algorithm filtering using a Kalman 
filter improves the results obtained, with standard deviations smaller than 0.5°: ci= 0.012°, 
c 2 = 0.227°, c 3 = 0.278° and c 4 = 0.189°. 

Future work will focus on exporting the know-how and the techniques developed in PAU-OR to the 
16 element PAU. 

Acknowledgements 

This work has been conducted as part of the award "Passive Advanced Unit (PAU): A Hybrid L- 
band Radiometer, GNSS-Reflectometer and IR-Radiometer for Passive Remote Sensing of the Ocean" 
made under the European Heads of Research Councils and European Science Foundation EURYI 
(European Young Investigator) Awards scheme in 2004, was supported by funds from the Participating 
Organizations of EURYI and the EC Sixth Framework Programme. It has also been supported by the 
Spanish National Research and EU FEDER Project TEC2005-06863-C02-01. 

References and Notes 

1. Font, J.; Lagerloef, G.S.E.; Le Vine, D.M.; Camps, A.; Zanife, O.Z. The determination of surface 
salinity with the European SMOS space mission. IEEE Transactions on Geoscience and Remote 
Sensing 2004, 42 (10), 2196-2205. 

2. Camps, A.; Font, J. ; Vall-llossera, M.; Gabarro, C; Corbella, I.; Duffo, N.; Torres, F.; Blanch, S.; 
Aguasca, A.; Villarino, E. The WISE 2000 and 2001 field experiments in support of the SMOS 
Mission: Sea Surface L-Band Brightness Temperature Observations And Their Application to 
Multi- Angular Salinity Retrieval. IEEE Transactions on Geoscience and Remote Sensing 2004, 42 
(4), 804-823. 

3. S0bjaerg, S.S.; Rotb0ll, J.M; Skou, N. The L-band Ocean Salinity Airborne Campaign. LOSAC 
ESA ISBN: SP-525, vol: LOSAC C, 2003, pp. 27-35. 

4. Reul, N. Contract ESA ITTl-4505/03/NL/Cb-Synergetic aspects and auxiliary data concepts for 
sea surface salinity measurements from space", technical report DOS, vol. 2004-03, IFREMER BP 
70, 29280 Plouzane, France. 2004. 



Sensors 2008, 8 



4411 



5. Sabia, R., Camps, A.; Vall-llossera, M.; Reul, N. Impact on Sea Surface Salinity Retrieval of Different 
Auxiliary Data Within the SMOS Mission. IEEE Transactions on Geoscience and Remote Sensing 
2006,44(11-1), 2769-2778. 

6. Camps, A. Telecommunications Engineering. "Passive Advanced Unit (PAU): A Hybrid L band 
Radiometer, GNSS-Reflectometer and IR-adiometer for Passive Remote Sensing of the Ocean", 
European Young Investigator Awards 2004, http://www.esf.org/activities/euryi/awards/2004.html 

7. Yueh, S.H. Modeling of Wind Direction Signals in Polarimetric Sea Surface Brightness 
Temperatures. IEEE Transactions on Geoscience and Remote Sensing 1997, 35, 6. 

8. Marchan-Hernandez, J.F.; Ramos-Perez, I.; Bosch-Lluis, X.; Camps, A. FPGA-based 
implementation of DDM-generator for GPS-reflectometry. Proceedings of the International 
Geoscience and Remote Sensing Symposium, USA 2006, 1180-1183. 

9. Marchan-Hernandez, J.F.; Ramos-Perez, I.; Bosch-Lluis, X.; Camps, A.; Rodriguez-Alvarez, N.; 
Albiol, D. PAU-GNSS/R, a Real-Time GPS-Reflectometer for Earth Observation Applications: 
Architecture Insights and Preliminary Results. Proceedings of the International Geoscience and 
Remote Sensing Symposium, Spain 2007. 

10. Marchan-Hernandez, J.F.; Rodriguez-Alvarez, N.; Camps, A. ; Bosch-Lluis, X.; Ramos-Perez, I.; 
Valencia, E. Correction of the Sea State Impact in the L-band Brightness Temperature by Means of 
Delay-Doppler Maps of Global Navigation Satellite Signals Reflected over the Sea Surface. IEEE 
Transactions on Geosceince and Remote Sensing 2008 (in press). 

11. Rodriguez-Alvarez, N., Camps, A.; Bosch-Lluis, X.; Ramos-Perez I.; Marchan-Hernandez, J.F. 
PAU-IR: Sea surface temperature retrieval using IR-radiometry and atmospheric modeling: 
Simulation and experimental results. International Geoscience and Remote Sensing Symposium, 
Spain 2007. 

12. Camps, A., Bosch-Lluis, X.; Ramos-Perez, I.; Marchan-Hernandez, J.F.; Izquierdo B.; Rodriguez- 
Alvarez, N. New Instrument Concepts for Ocean Sensing: Analysis of the PAU-Radiometer. IEEE 
Transactions on Geosceince and Remote Sensing 2007, 45 (10), 3180-3192. 

13. Ramos-Perez, I., Bosch-Lluis, X.; Camps, A.; Marchan-Hernandez, J.F.; Prehn, R. Design of a 
compact dual-polarization receiver for pseudo-correlation radiometers at L-band. Proceedings of 
the International Geoscience and Remote Sensing Symposium, USA, 2006, 1172-1175. 

14. Bosch-Lluis, X.; Camps, A.; Marchan-Hernandez, J. F.; Ramos-Perez, I.; Prehn, R. FPGA-based 
implementation of a polarimetric radiometer with digital beamforming. Proceedings of the 
International Geoscience and Remote Sensing Symposium, USA 2006, 1 176-1 179. 

15. Kalman, R.E. A new approach to linear filtering and prediction problems. Transactions of the 
ASME-Journal of Basic engineering 1960, 82, 35-40. 

16. Smith, N.; Lefebvre, M. The global ocean data assimilation experiment (GODAE). in Proc. Int. 
Symp. Monitoring Oceans in the 2000s: An Integrated Approach, France 1997. 

17. Camps, A., Rodriguez-Alvarez, N.; Bosch-Lluis, X.;. Marchan, J.F; Ramos Perez, I.; Segarra, M.; 
Sagues, LI.; Tarrago, D.; Cunado O. PAU in SeoSAT: A Proposed Hybrid L-band Microwave 
Radiometer/GPS Reflectometer to Improve Sea Surface Salinity Estimates from Space. 10 th 
Specialist Meeting on Microwave Radiometry and Remote Sensing of the Environment, Italy 2008. 



Sensors 2008, 8 



4412 



18. Van Vlecr, J.H.; Middleton, D. The Spectrum of Clipped Noise. Proceedings of the IEEE 1966, 54 
(1), 2-19. 

19. Bosch-Lluis, X.; Ramos-Perez, I.; Camps, A.; Marchan-Hernandez, J.F.; Rodriguez-Alvarez, N. 
Valencia, E.; Vericat, M.; Guerrero, MA. The Impact of the Number of Bits in Digital 
Beamforming Real Aperture and Synthetic Aperture Radiometers. 10 th Specialist Meeting on 
Microwave Radiometry and Remote Sensing of the Environment, Italy 2008. 

20. De Paula, A.; Pieper, R.J. A more complete analysis for subnyquist band-pass sampling. System 
Theory, 1992. Proceedings of the SSST/CSA 92, 1992, 20-24. 

21. Yueh, S.H. Estimates of Faraday rotation with passive microwave polarimetry for microwave 
remote sensing of Earth surfaces. IEEE Transactions on Geosceince and Remote Sensing 2000, 38 
(5), 2434-2438. 

22. Camps, A., Aguasca, A.; Bosch-Lluis, X.; Marchan-Hernandez, J.F.; Ramos-Perez, I.; Rodriguez- 
Alvarez, N.; Bou, F.; Ibanez, C. Measurements and First Results of the One-Receiver PAU 
Instruments. International Geoscience and Remote Sensing Symposium, Spain 2007. 



© 2008 by the authors; licensee Molecular Diversity Preservation International, Basel, Switzerland. 
This article is an open-access article distributed under the terms and conditions of the Creative 
Commons Attribution license (http://creativecommons.Org/licenses/by/3.0/). 



